Storage of energy in confined quantum systems 
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Using the non-perturbative method of dressed states intro- 
duced in previous publications, we study the evolution of a 
confined quantum mechanical system embedded in a ohmic 
environment. Our approach furnishes a theoretical mecha- 
nism to control inhibition of the decay of excited quantum 
systems in cavities, in both weak and strong coupling regimes. 
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In recent publications [1,2] a method [dressed states) 
has been introduced that allows a non-perturbative ap- 
proach to situations where a material body (an atom or a 
Brownian particle) is coupled to a field, provided that the 
interaction between the parts of the system can be ap- 
proximated by a linear coupling. From now on we use the 
term "particle" in a general sense, refering to an atom, a 
Brownian particle a molecule, etc.... Our dressed states 
can be viewed as a rigorous version of the semiqualita- 
tive idea of dressed atom introduced originally in ref . [3] , 
which can be constructed in reason of the linear charac- 
ter of our problem. More precisely, the method applies 
for all systems that can be described by an Hamiltonian 
of the form, 
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where the subscript refers to the "material body" and 
= 1,2, ...N refer to the harmonic environment modes. 
A Hamiltonian of this type, describing a linear coupling 
of a particle with an environment, has been used in [4] to 
study the quantum Brownian motion of a particle with 
the path-integral formalism. The limit ^ oo in Eq.(l) 
is understood. We investigate also the behaviour of the 
system as a function of the strenght of the coupling be- 
tween the particle and the bath. In particular we give 
explicitly non-perturbative formulas for the decay prob- 
ability of the particle in the weak and strong coupling 
regimes. 

We consider the problem of a particle qo in the har- 
monic approximation having (bare) frequency luq coupled 
to A^ other oscillators qi of frequencies uji, i — 1,2, ...N. 
In the limit N —^ oo we recover our situation of the 
coupling particle-bath (field) after redefinition of diver- 
gent quantities. The bilinear Hamiltonian (1) can be 
turned to principal axis by means of a point transfor- 
mation, = V^Qr Pf^ = t^^Pr; IJ. = (0, {fc}), k = 
1,2,...,A^; r = 0, ...A, performed by an orthonormal 
matrix T = (i^. The subscript ^ = Q and fi = k refer re- 
spectively to the particle and the harmonic modes of the 



bath and r refers to the normal modes. In terms of nor- 
mal momenta and coordinates, the transformed Hamil- 
tonian in principal axis reads, H = j ^r^oi-^r + ^rQf): 
where the fi^'s are the normal frequencies corresponding 
to the possible collective stable oscillation modes of the 
coupled system. The matrix elements are given by [1] 
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with the condition, 
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We take Ck = rj{ujk)^ . In this case the environment is 
classified according io n > l,n = l, orn< 1, respec- 
tively as supraohmic, ohmic or subohmic. For a subohmic 
environment the sum in Eq.(3) is convergent and the fre- 
quency ujQ is well defined. For ohmic and supraohmic 
environments the sum in the right hand side of Eq.(3) di- 
verges what makes the equation meaningless as it stands, 
a renormalization procedure being needed. In this case, 
by a straightforward generalization of the method de- 
scribed in [1] we can define a renormalized frequency uj, 
by means of a counterterm Suj'^, 
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in terms of which Eq.(3) becomes. 
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with the notation E[n\ standing for the smallest integer 
containing n. We see that in the limit N oo the above 
procedure is exactly the analogous of naive mass renor- 
malization in Quantum Field Theory: the addition of a 
counterterm — Jcj^gg allows to compensate the infinuty of 
ujq in such a way as to leave a finite, physically meaninful 
renormalized frequency u. This simple renormalization 
scheme has been originally introduced in ref. [5]. 

To proceed, we take the constant rj as rj = ^/2gAu!, 
Auj being the interval between two neighbouring bath 
frequencies (supposed uniform) and where g is some con- 
stant [with dimension of {frequency)^~^]. For reasons 
that will become apparent later, we restrict ouselves to 



the physical situations in which the whole system is con- 
fined to a cavity of diameter L, and the environment 
frequencies Uk can be writen in the form 



Wk = ^kir/L, fc = l,2, ... . 

Then using the formula, 
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and restricting ourselves to an ohmic environment, Eq.(5) 
can be written in closed form. 
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The solutions of Eq.(8) with respect to Q give the spec- 
trum of eigcnfrequencies corresponding to the collec- 
tive normal modes. 

The transformation matrix elements turning the material 
body-bath system to principal axis is obtained in terms 
of the physically meaningful quantities fi,.. lj, after some 
rather long but straightforward manipulations analogous 
as it has been done in [1]. They read, 
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To study the time evolution of the particle, we start 
from the eigenstates of our system, represented by the 
normalized eigenfunctions. 
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where stands for the n^-th Hermite polynomial and 
To is the normalized vacuum eigenfunction. We intro- 
duce (ire.s.sed coordinates q'^ and {q'^} for, respectively the 
dressed atom and the dressed field, defined by [1] , 
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valid for arbitrary L and where cj^ = {tj, cjj}, and let us 
define for a fixed instant the complete orthonormal set of 
functions [1], 
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where = q^, q^, uj^ = {w, Wj}. Note that the groimd 
state To in the above equation is the same as in Eq.(lO). 



The invariance of the ground state is due to our definition 
of dressed coordinates given by Eq. (11). Each function 
V'kqki- - describes a state in which the dressed oscil- 
lator q'^ is in its — th excited state. Using Eq.(ll) 
the functions (12) can be expressed in terms of the nor- 
mal coordinates Qr- But since (10) is a complete set of 
orthonormal functions, the functions (12) may be writ- 
ten as linear combinations of the eigenfunctions of the 
coupled system, and conversely. 

Let us call T\ = V'ioo...o(m)o... the configuration in which 
only the dressed oscillator (the dressed particle) is in 
the first excited level. Then it is shown in [1] the following 
expression for the time evolution of the first-level excited 
dressed particle <7q, 
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Prom Eq.(13) we see that the initially excited dressed 
particle naturally distributes its energy among itself and 
all other dressed oscillators (the environment) as time 
goes on, with probability amplitudes given by the quan- 
tities fo„{t) in Eq.(13). In Eq.(13) the coefficients f''{t) 
have a simple interpretation: /°°(t) and /°'(t) are re- 
spectively the probability amplitudes that at time t the 
dressed particle still be excited or have radiated a quan- 
tum of frequency hwi. We see that this formalism allows 
a quite natural description of the radiation process as a 
simple exact time evolution of the system. In the case of 
a very large cavity (free space) our method reproduces for 
weak coupling the well-known perturbative results [1,2]. 
Nevertheless it should be remarked that due to our ap- 
proximations to be done below for the confined system, 
we can not recover free space results from our confined 
system. 

Next we study the time evolution process of the con- 
fined dressed particle when it is prepared in such a way 
that initially it is in its first excited state. We shall con- 
sider the situation of the particle confined in a cavity of 
diameter L. 

Let us now consider the ohmic system in which the 
particle is placed in the center of a cavity of diameter 
L, in the case of a small L, i.e. that satisfies the con- 
dition of being much smaller than the coherence lenght, 
L « 2c/ g. We note that from a physical point of view, 
L stands for either the diameter of a spherical cavity or 
the spacing between infinite paralell mirrors. To obtain 
the eigcnfrequencies spectrum, we remark that from a 
graphical analysis of Eq.(8) it can be seen that in the 
case of small values of L, its solutions are very near 
the frequency values corresponding to the asymptots of 
the curve cot(^), which correspond to the environment 
modes uji ~ i2TTc/L, except from the smallest cigcnfre- 
quency VIq. As we take larger and larger solutions, they 
are nearer and nearer to the values corresponding to the 



asymptots. For instance, for a value of L of the order 
of 2 X 10~^m and a) ~ 10^°/s, only the lowest eigen- 
frequency Oq is significantly different from the field fre- 
quency corresponding to the first asymptot, all the other 
eigenfrequencies 17^, k — 1,2,... being very close to the 
field modes fc27rc/L. For higher values of uj (and lower 
values of L) the differences between the eigenfrequencies 
and the field modes frequencies are still smaller. Thus 
to solve Eq.(8) for the larger eigenfrequencies we expand 
the function cot(^) around the values corresponding to 
the asymptots. We write. 
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fife = — (A; + efc), /c = l,2, . 
with < Cfc < 1, satisfying the equation, 
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But since for a small value of L every Ck is much smaller 
than 1, Eq.(15) can be linearized in e^, giving, 
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Eqs.(14) and (16) give approximate solutions to the 
eigenfrequencies fife, fc=l,2.... 

To solve Eq.(8) with respect to the lowest eigenfre- 
quency flo, let us assume that it satisfies the condition 
iloL/2c « 1 (we will see below that this condition is 
compatible with the condition of a small L as defined 
above). Inserting the condition QqL/2c << 1 in Eq.(8) 
and keeping up to quadratic terms in we obtain the 
solution for the lowest eigenfrcquency fio, 
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Consistency be- 
tween Eq.(17) and the condition fioL/2c << 1 gives a 
condition on L, 



(18) 

Let us consider the situations of weak coupling {g <C w) 
and of strong coupling {g lo), and let us consider the 
situation where the dressed material body is initially in 
its first excited level. Then from Eq.(13) we obtain the 
probability that it will still be excited after a ellapsed 
time t, 
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a) Weak eoupling 
In the case of weak coupling a physically interesting sit- 
uation is when interactions of electromagnetic type are 
involved. In this case, we take g = au, where a is the 
fine structure constant, a — 1/137. Then the factor / 
multiplying 2c/ g in Eq.(18) is ~ 0.07 and the condi- 
tion L <C 2clg is replaced by a more restrictive one, 
L ^ 0.07(2c/,g). For a typical infrared frequency, for 
instance w ^ 2, x 10^^/s, our calculations are valid for 
a value of X, L <C 10~^m. 

From Eqs.(9) and using the above expressions for the 
eigenfrequencies for small L, we obtain the matrix ele- 
ments, 
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To obtain the above equations we have neglected the 
corrective term e^, from the expressions for the eigen- 
frequencies fife. Nevertheless, corrections in Cfe should be 
included in the expressions for the matrix elements t^., in 
order to avoid spurious singularities due to our approxi- 
mation. 

Using Eqs.(20) in Eq.(19), we obtain 
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where we have introduced the dimensionless parameter 
S = Lg/2c, corresponding to a small value of L and we re- 
member that the eigenfrequencies are given by Eqs.(14). 
(16) and (17). As time goes on, the probability that the 
material body be excited attains periodically a minimum 
value which has a lower bound given by. 
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For instance, for a frequency lj of the order uj ^ 4.00 x 
W^/s (in the red visible), and L ~ 1.0 x 10" TO, remem- 
bering our weak coupling anzatz g = ua, we see from 
Eq.(22) that the probability that the material body be 
at any time excited will never fall below a value ^ 0.97, 
or a decay probability that is never higher that a value 
^ 0.03. It is interesting to compare this result with ex- 
perimental observations in [6] , where stability is found for 
atoms emitting in the visible range placed between two 
parallel mirrors a distance L = 1.1 x 10~^to apart from 
one another [6] . For cavities small enough the probability 
that the material body remain excited as time goes on, 
oscillates with time between a maximum and a minimum 
values and never departs significantly from the situation 
of stability in the excited state. 



Nevertheless, in the weak couphng regime for any fixed 
value of the emission frequency, the probability that 
the particle remain excited can be considerably changed. 
This is a general phenomenon. We see from Eq. (22) 
that mathematically there are two values of the cavity 
size, corresponding to 5 = and 5 = 15/147r which corre- 
spond to the probability one of the particle be excited at 
any time. Also, we can see from Eq.(22) that the minimal 
value of Min(|fO°(t)p) is attained for 5 = 15/287r. The 
size of the cavity corresponding to this value of S is ex- 
actly half the (non-zero) size of the cavity allowing from 
a mathematical point of view, permanent excitation of 
the particle. However we should not give physical mean- 
ing to the portion of the parabole (22) from S = 15/287r 
to infinity, it will correspond for values S > 15/147r to 
probabilities larger than one and is physically meaning- 
less. The physically meaningful part of the parabole (22) 
corresponds to < 5 < 15/287r. Our conclusion is that 
varying S from to J = 15/287r we attain at the mid- 
dle of the curve the point of maximum emission, where 
Min(|f°'^(t)|^) w 0.55. In the weak coupling regime, this 
is the minimal possible probability to the particle remain 
excited inside the cavity for any emmission frequency. 

b) Strong coupling 
In this case it can be seen from Eqs (16), (17) and (18) 
that flo Pi ui and we obtain from Eq.(9), 

Using Eqs. (23) in Eq.(19), we obtain for the probability 
that the excited system be still at the first energy level 
at time t, the expression, 

^ fe=i 
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+ ^^' E ^oos{n,-ni)t. (24) 

We see from (24) that as the system evolves in time, the 
probability that the material body be excited, attains 
periodically a minimum value which has a lower bound 
given by. 

The behaviour of the system in the strong coupling 
regime is completely different from weak coupling be- 
haviour. The condition of positivity of (25) imposes for 
fixed values of g and ui an upper bound for the quantity 
S, 5maxi which corresponds to an upper bound to the 
diameter L of the cavity, L^^ax (remember 5 = Lg/2c). 
Values of S larger than S^nax, or equivalently, values of 
L larger than L^^x are unphysical (correspond to nega- 
tive probabilities) and should not be considered. These 



upper bounds are obtained from the solution with re- 
spcict to (5, of the inequality, Min((5) > 0. We have 
Min((5) > or Min((5) = 0, for respectively S < S^ax or 
S — Smax- The solution of the equation Min((5) = gives, 

Smax = + « 0.37. For d = 6max, the min- 

imal probability that the atom remain stable vanishes. 
We see comparing with the results of the previous subsec- 
tion, that the behaviour of the system for strong coupling 
is rather different from its behaviour in the weak coupling 
regime. For appropriate cavity dimensions, which are of 
the same order of those ensuring stability in the weak 
coupling regime, we ensure for strong coupling the com- 
plete decay of the particle to the ground state in a small 
ellapsed time. 

We have presented in this note a non-perturbative 
treatment of a confined ohmic quantum system consist- 
ing of particle (in the larger sense of a " material body" , 
an atom or a Brownian particle) coupled to an environ- 
ment modeled by non- interacting oscillators. One possi- 
ble conclusion is that by changing conveniently the phys- 
ical and geometric parameters (the emission frequency, 
the strenghness of the coupling and the size of the con- 
fining cavity) our formalism theoretically allows a control 
on the rate of emission and of the energy storage capac- 
ity (perhaps information?) in the cavity. Depending on 
the strenghness of the coupling, the emission probability 
ranges from a complete stability to a very rapid decay. 
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